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Abstract. The gamma-ray observations of molecular clouds associated with supernova remnants are considered 
one of the most promising ways to search for a solution of the problem of cosmic ray origin. Here we briefly 
review the status of the field, with particular emphasis on the theoretical and phenomenological aspects of the 
problem. 


1 Introduction 

In a pioneering paper. Black and Fazio (Tj proposed to 
measure the mass of Molecular Clouds (MCs) from their 
gamma-ray emission. This implicitly assumes the exis¬ 
tence of an ubiquitous and uniform sea of Cosmic Rays 
(CRs) that pervades the whole Galaxy. Under this assump¬ 
tion the intensity of CRs in any given MC is known and 
its gamma-ray emission, dominated by the contribution 
from the 7r°-decay channel, simply scales as F y oc 2 , 
where M c i is the mass of the cloud and d its distance (see 
e.g. 0). Thus, if the distance is also known, the mass of 
the cloud can be estimated by measuring F y . 

In fact, the mass of a MC can be estimated also from 
the intensity of the CO line emission 0, and thus the ar¬ 
gument above can be reversed and one can then infer, from 
gamma-ray observations, the intensity of CRs into the MC 
under examination (see e.g. 00). For this reasons MCs 
have been sometimes referred to as cosmic ray barome¬ 
ters, because they can serve as probes to estimate the CR 
intensity (and thus their pressure) at various locations in 
the Galaxy 0. 

Both the aforementioned approaches have limitations. 
The assumption of a rough spatial homogeneity of CRs is 
justified only on large spatial scales, but is most likely vi¬ 
olated in the vicinity of particle accelerators, due to the 
injection of freshly accelerated CRs in the surrounding in¬ 
terstellar medium 00- On the other hand, mass esti¬ 
mates based on the observations of the CO line emission 
rely on the poorly determined CO-toTU conversion fac¬ 
tor, Aco, and this adds an uncertainty in the use of MC as 
CR barometers lUOil . However, excesses above the CR sea 
of one or two order of magnitudes can be easily realized 
in the proximity of CR sources, clearly overwhelming the 
uncertainty on Aco (see Table 1 in 0). To conclude, MCs 
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can be definitely used as probes to detect the presence of 
large and localized enhancements in the CR intensity, and 
thus the presence of CR sources. Clearly, this fact has cru¬ 
cial implications in the quest for the origin of Galactic CRs 
(see for an extended review). 

Excesses in the CR intensity have been indeed detected 
by means of gamma-ray observations of MCs (e.g. US¬ 
ED), and in most cases the MCs are located in the vicinity 
of, or partially overlap with SuperNova Remnants (SNRs), 
the putative sources of CRs. Remarkably, for two SNRs 
(W28 and W44) part of the gamma-ray emission comes 
from a region located clearly outside of the SNR shell. In 
this case, the observations can be interpreted as the result 
of the escape of CRs from the SNR shock, their propaga¬ 
tion up to the cloud, and their interaction with the dense 
gas na. In the following Section we will describe how 
the interpretation of these observations might lead to an 
indirect proof for the fact that SNRs indeed accelerate CR 
hadrons, and thus reinforce the belief that they might be 
the sources of Galactic CRs. Also, a constrain on the prop¬ 
agation properties of CRs in the interstellar magnetic field 
can be obtained from such observations. 

An implicit assumption done so far is that CRs pene¬ 
trate unimpeded MCs, i.e., the CR spectrum immediately 
outside of the MC is identical to that inside it. In fact, this 
might not be the case, due to the interplay of CR spatial 
transport at the MC bondary and energy losses suffered by 
CRs inside MCs fl8lH9l . We will briefly discuss this issue 
in Section 3. Finally, in Section 4 we will present a multi¬ 
wavelength approach aimed at extracting the CR spectrum 
from MCs over more than six orders of magnitude in en¬ 
ergy: from the sub-MeV to multi-TeV domain. 

2 Runaway cosmic rays 

Though the details of the escape of CRs from SNR shocks 
are still debated (see e.g. j20l . or EQ and references 
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Figure 1 . Fit of the gamma-ray emission from three molecu¬ 
lar clouds located in the vicinity of the supernova remnant W28. 
Isotropic diffusion of cosmic rays has been assumed. Data points 
refer to HESS observations lH3l . Figure from |j22l . 


therein), it is generally believed that the highest energy 
particles (we focus here on protons) are released first, and 
particles of lower energy are released gradually as the 
shock speed decreases. Once released, relativistic pro¬ 
tons produce gamma-rays due to their inelastic interac¬ 
tions with the surrounding interstellar medium 

This is of particular interest if a massive MC is located 
in the vicinity of the CR source, given that the amount of 
gamma rays produced in CR interactions scales with the 
mass of the target material. In this context, the gamma-ray 
emission produced by a MC illuminated by CRs escaped 
from a nearby SNR has been computed in (9j QJI for an 
isotropic CR diffusion coefficient. These results have been 
then applied to the massive MCs located in the vicinity of 
the SNR W28 (22) and a good fit to the observed gamma 
ray emission has been presented (see Fig. [T|i. 

The assumption of isotropy of the CR diffusion coef¬ 
ficient D implies that within a diffusion length from the 
SNR, R,/ ~ V£> t esc , the CR intensity is roughly /cr oc 
r/Es^/Rj, where £sn ~ 10 51 erg is the supernova explo¬ 
sion energy and 77 is an efficiency that represents the frac¬ 
tion of such energy converted into CRs. If SNRs are the 
sources of CRs it must be 77 « 10%. The time since CRs 
left the remnant is t esc which, for CRs of sufficiently large 
energy and for sufficiently old SNRs, is comparable to the 
SNR age: t ese ~ t age G7). For the case of W28 t age is of the 
order several 10 4 years. The gamma-ray flux from MC of 
known mass M c \ is then F y oc jcRM c i/d 2 where d ~ 2 kpc 
is the distance to W28. If the diffusion coefficient is ex¬ 
pressed in units of the typical galactic one D = xD ga i, 
with D ga i 10 28 (£710 GeV ) 0 5 cm 2 /s, one can see that the 
gamma-ray flux scales proportionally to the ratio tj/x 3 ^ 2 . 


Thus, with this kind of studies one can not only indirectly 
reveal the acceleration of hadronic CRs at SNR shocks, 
and possibly prove (or falsify) the SNR hypothesis for the 
origin of CRs, but also constrain the CR diffusion coeffi¬ 
cient in the vicinity of the SNR. The fit in Fig. 1 has been 
obtained assuming r//x 2 ^ 2 ~ 20 , which implies that for any 
plausible value of the CR acceleration efficiency 77 < 1, the 
CR diffusion coefficient in the vicinity of the SNR has to 
be significantly smaller than the typical galactic one (i.e. 
X <s 1) ll22l . Similar conclusions have been independently 
obtained also in (23W25) . 

However, the constrains on the diffusion coefficient 
change significantly if the assumption of isotropy of the 
CR diffusion is relaxed. In fact, over scales smaller than 
the galactic magnetic field correlation length, the diffusion 
of CRs is likely to be more effective along field lines rather 
than perpendicular to them (see e.g. lf26M28il ). This is 
shown in Fig. 2 (from (26l ). where the expected CR over¬ 
density above the galactic CR sea is mapped for a region 
surrounding a SNR. The supernova explosion is assumed 
to have happened in the centre of the plot, and an isotropic 
(left panel) and anisotropic (right panel) diffusion coeffi¬ 
cient has been adopted. For an anisotropic diffusion coef¬ 
ficient runaway CRs remain focussed along the field lines, 
and thus their intensity there is correspondingly enhanced. 

In Ii26l the gamma-ray emission from the MCs in the 
W28 region has been fitted by assuming a fully anisotropic 
diffusion of CRs. To fit data, a much larger diffusion co¬ 
efficient has to be assumed with respect to the case of 
isotropic diffusion, namely D|| ~ 10 2,s (£/10 GeV ) 05 cm 2 /s 
(which naively would correspond to^f ~ 1 , though a direct 
comparison between an isotropic and a parallel diffusion 
coefficient would require a more sophisticated discussion). 

Summarising, the interpretation of the gamma-ray ob¬ 
servations of MCs illuminated by CRs accelerated at 
nearby SNRs dramatically depends on the assumptions 
made on the properties of CR diffusion. Both observations 
with future facilities such as the Cherenkov Telescope Ar¬ 
ray (29) and thorough theoretical investigations (30) are 
needed in order to solve this issue. 

3 Penetration of cosmic rays in molecular 
clouds 

An assumption implicitly done so far in this paper is that of 
unimpeded penetration of CRs into MCs, i.e., the spectrum 
and intensity of CRs inside and immediately outside of the 
cloud are identical. However, there are reasons to believe 
that this might not be the case ITEl [T9! 1321434) . If so, the 
gamma-ray emission from MCs might be correspondingly 
affected USES 

The main difficulty of studying from a theoretical point 
of view the problem of CR penetration in MCs resides in 
its strong non-linearity: CRs themselves excite the mag¬ 
netic turbulence that in turn scatters them and obstructs 
their penetration into the cloud. The plasma instability re¬ 
sponsible for the generation of the turbulence is the cosmic 
ray streaming instability (36). Such instability operates at 
the border of the cloud only, at the transition between ion- 
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Figure 2. Cosmic ray over-density around a typical supernova remnant (see |26| for details) for a particle energy of E = 1 TeV 
at a time t = 10 kyr after the explosion. The left panel refers to an isotropic diffusion coefficient of cosmic rays equal to D = 
5 x 10 26 (£710 GeV) 05 cm 2 /s, while the right panel refers to an anisotropic diffusion scenario with Dy = 10 28 (£/10GeV)° 5 cm 2 /s and a 
magnetic field coherence length equal to 20 pc. The black cross marks the a position at which the CR over-density is equal in the two 
panels. 


ized and neutral gas. The reason for that is that in the neu¬ 
tral medium inside the cloud the turbulence is immediately 
damped by ion-neutral friction, making the streaming in¬ 
stability ineffective li37ll . 

However, it was recently pointed out in m that the 
solution of the problem has a universal character and does 
not depend on the details of the plasma instability that de¬ 
termines the properties of CR propagation. To understand 
this, let us recall that the intensity of low energy CRs in¬ 
side the cloud is reduced because particles undergo severe 
ionization energy losses there l38l . Therefore, to match 
the interstellar spectrum of CRs far away from the cloud, a 
spatial gradient in the CR intensity must appear at the out¬ 
skirts of the cloud. This is the region where CR streaming 
instability operates, and where CRs generate the magnetic 
turbulence that scatters them and determines their trans¬ 
port properties. The crucial result obtained in lfl9l is that 
such a gradient that forms in the outskirts of the cloud 
acts as a self-regulating buffer for the flux of CRs enter¬ 
ing the cloud. If streaming instability is very efficient (or 
inefficient) then the extension of the zone characterized by 
the gradient changes accordingly, shrinking (or expand¬ 
ing) and keeping in this way the flux of the CRs into the 
cloud constant and equal to the universal value v A fo, where 
v A is the Alfven speed in the diffuse interstellar medium 
(the velocity of the magnetic perturbations of the field) 
and /o the interstellar spectrum of CRs. Thus, the spec¬ 
trum of CRs inside the MC can be described in terms of a 
characteristic energy 03: 

(i » 


where Nh is the cloud gas column density and a ~ 2.58 is 
the slope of the ionization energy loss time as a function 
of the particle momentum: r IO „ oc p". Particles of energy 
larger than Eb r are unaffected by energy losses, thus the 
CR spectrum inside and outside of the MC is identical in 
this energy domain. On the other hand, particles with ener¬ 
gies smaller than Et,,- suffer severe ionization losses so that 
the CR spectrum into the cloud is affected significantly. If 
the CR spectrum in the interstellar medium is a power law 
in momentum /o oc p ', then the spectrum inside the cloud 
at particle energies smaller than E hr reads: 

f c oc p" -3 s < a - 3 

fc cx p- s s>a- 3 (2) 

Remarkably, for a CR spectrum in the interstellar medium 
fo(p) ~ p~ s with s = a-3 ~ 0.42 the slope of the spectrum 
of CR is identical inside and outside of the cloud. CR 
spectra inside and outside of a MC of column density 3 x 
10 21 cm~ 2 are plotted in Fig. 3. 

To conclude, we note that a break in the spectrum in 
the * 100 MeV range (see expression for Eq. 1) range 
would not affect significantly the gamma-ray luminosity 
of the cloud, because the threshold for neutral pion pro¬ 
duction in proton-proton interactions is at a larger energy, 
namely ^ 280 MeV. A possibly unique exception to this is 
the Sgr B2 MC complex, the most massive MC complex 
in the Galaxy, characterised by an extremely large mass 
(w 10 7 M q ) and gas column density (up to 2.5 x 10 24 cm~ 2 
1401 ). The gamma-ray emission from such cloud reveals 
a break in the CR parent spectrum at a particle energy of 
~ 4 GeV lf39l . which might be a manifestation of Eq. 1. 
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Figure 3. Spectra of cosmic rays inside a molecular cloud of column density Nh = 3 x 10 21 cnT 2 (solid lines) assuming that the spectra 
far for the cloud are given by power law in momentum of slope s (dashed lines). The dot-dashed line shows the eigenfunction of the 
problem. 


4 A multi-wavelength approach: from 
millimeter to TeV observations 

The importance of studying the amount of penetration of 
low energy CRs into MCs is connected to the fact that 
CRs are the only ionizing agents capable of penetrating 
large gas column densities and thus regulate the ionization 
fraction of dense MC cores. The ionization fraction in 
turns determines the level of coupling between the mag¬ 
netic field and the gas, thus playing a pivotal role in the 
dynamical evolution of MCs EE). 

After an early interest on this issue (see e.g. osicsa- 
(34) , theoretical investigations of the level of penetration 
of low energy cosmic rays into molecular clouds became 
quite rare, mainly due to the difficulties connected to the 
measurements of both the ionization level of clouds (but 
see e.g. ED for a notable exception) and of the interstel¬ 
lar cosmic ray spectrum at energies low enough (~ MeV) 
to be relevant for gas ionization (e.g. Ii42l ). The break¬ 
through came during the last decade when observations 
of the infrared H.J line and of the DCO + /HCO + millime¬ 
ter lines ratio from molecular clouds provided us with the 
measurement of the cosmic ray ionization rate from a quite 
significant sample of clouds (e.g. (43), for reviews see 
mmi This revived the interest in theoretical studies of 
the cosmic ray penetration into clouds lf38l |46l l. In addi¬ 
tion to that, the Voyager 1 probe, having reached a distance 
of 122 astronomical units from the Earth, reported on the 
measurement of the interstellar spectrum of CRs, i.e. the 
spectrum unaffected by solar modulation, down to ener¬ 
gies of several MeV per nucleon (Krimigis et al. 2013, 
Stone et al. 2013). These measurements reduce dramati¬ 
cally the uncertainty on our knowledge of the CR spectrum 


in the interstellar medium (at least the local one), which 
however remains unconstrained below energies of several 
MeV, that are still quite relevant for ionization (38). 

Also in this context, associations of SNRs with MCs 
are of prime interest. If a MC interacts with a SNR shock, 
one naturally expects to observe, besides an enhanced 
(hadronic) gamma-ray emission, also an increased CR ion¬ 
isation rate. The relative enhancements would depend on 
the shape of the spectrum of CRs accelerated at the SNR 
shock itself, since * MeV CRs are responsible for ioni¬ 
sation, while ~ GeV to TeV CRs are responsible for the 
production of gamma-rays. 

Over the last several years, CR ionisation rates at the 
level of > 10 15 s“ 1 have been reported for MCs interacting 
with the shocks of the SNRs IC 443, W51c, and W28 (see 
mm, respectively). Since these values are well above 
the ionisation rate one would expect in an isolated cloud 
immersed in the galactic CR sea, such observations might 
constitute the smoking gun for the acceleration of low en¬ 
ergy CRs at SNR shocks. However, the relative contribu¬ 
tion of CR protons and electrons to the ionisation rate still 
remains undetermined 1351 . Also the contribution to ion¬ 
isation from X-ray photons from the SNR interior needs 
to be properly assessed (unlike UV photons. X-rays can 
penetrate quite large column densities of gas). 

If the origin of the enhanced ionisation revealed in the 
above mentioned clouds will be found to have an hadronic 
origin, then a combination of low and high energy obser¬ 
vations will allow to constrain the spectrum of CRs freshly 
injected in the interstellar medium by a SNR shock over an 
energy interval of unprecedented breadth: from the MeV 
to the multi-TeV domain. 
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As an example, let us now consider the MC in inter¬ 
action with the shock of the SNR W28. A ionisation rate 
of > 2 x 10~ 15 s _1 has been measured there iBTI . Also, 
gamma-rays of very likely hadronic origin have been de¬ 
tected by FERMI and H.E.S.S. in the GeV and TeV do- 
mian, respectively ED IH. Then, let us fit the gamma-ray 
emission by assuming the underlying CR proton spectrum 
to be a power law in momentum. This requires a spectral 
slope of 4.8, or 4np 2 N(p) oc p~ 2 S . The normalisation of 
N(p) can be derived from the measurement of the gamma- 
ray flux of the MC and from the values of the MC mass 
and distance, which have been estimated to be equal to 
5 x 1O 4 M 0 and 2 kpc, respectively. Once the normalisa¬ 
tion is obtained, the expected ionisation rate due to CR 
protons can be computecQ This of course will depend on 
the extension of the CR spectrum to low energies, which 
are more relevant for ionisation. If the CR proton spectrum 
is abruptly truncated below a particle energy of 280, 100, 
and 60 MeV, the corresponding ionisation rates would be 
»4x 10~ 16 , 10^ 15 , and 2 x 10~ 15 s _I , respectively. This 
implies that in order to be consistent with the measured 
ionisation rate reported in ED, the CR proton spectrum 
cannot continue with the same slope below particle en¬ 
ergies of « 50 - 100 MeV. A break or a cutoff must be 
present in the CR spectrum. Also, it seems that there is 
not much room left for electron and/or X-ray induced ion¬ 
isation of the gas, unless the CR proton spectrum is cutted 
off abruptly at a particle energy just below or of the order 
of the threshold for neutral pion production (w 280 MeV). 

These simple estimates exemplify how constrains on 
the shape of the CR spectrum can be obtained by combin¬ 
ing low and high energy observations of MCs interacting 
with SNRs. This would shed new light on the process of 
particle acceleration at SNR shocks, and on the transport 
of CRs after their injection into the interstellar medium. A 
complete understanding of these aspects of CR physics is 
needed in order to verify (or not) the SNR hypothesis for 
the origin of galactic CRs. 
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